Multiple condensed phases in attractively interacting Bose systems 
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We investigate a Bose gas with finite-range interaction using a scheme to eliminate unphysical 
processes in the T-matrix approximation. In this way the corrected T-matrix becomes suitable to 
calculate properties below the critical temperature. For attractive interaction, an Evans-Rashid 
transition occurs between a quasi-ideal Bose gas and a BCS-like phase with a gaped dispersion. The 
gap decreases with increasing density and vanishes at a critical density where the single-particle dis- 
persion becomes linear for small momenta indicating Bose-Einstein condensation. The investigation 
of the pressure shows however, that the mentioned quantum phase transitions might be inaccessible 
due to a preceding first order transition. 
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The multiple-scattering correlations collected in terms 
of the T-matrix are a powerful tool to describe properties 
of interacting many-body systems. Compared to the T- 
matrix of the scattering theory, the many-body T-matrix 
contains pairs of medium-dependent intermediate prop- 
agators. They describe the influence of the surround- 
ing medium and their selfconsistent determination is a 
challenging task. It was first studied in dense Fermi sys- 
tems, where taking the particle-particle correlations as 
well as hole-hole correlations into account a combination 
of self-consistent propagators appear. The application 
of this fully self-consistent approach called the Galitskii- 
Feynman approximation below the critical temperature 
of condensation is questionable. Considering the case of 
cooper pairing the break down of this approximation can 
be observed by the fact that although the divergence of 
the T-matrix leads to the correct critical temperature 
the Galitskii-Feynman approximation does not lead to 
the correct gap dispersion for the condensed phase 

It was observed by Kadanoff and Martin Q and used 
later on [H, H, 0, [H that an asymmetric breaking of the 
selfconsistency in the T-matrix, such that one propagator 
is used selfconsistently and the other non-selfconsistently, 
leads immediately to the correct gap dispersion while the 
symmetric selfconsistent one does not. This so called 
Prange paradox [|| remains puzzling since a seemingly 
worse approximation leads to better results. Recently it 
was found that the reason for that could be unphysical 
repeated collisions [l|. 

While in the thermodynamic limit the contributions 
of these or similar unphysical processes usually vanish, 
since each channel contributes separately with the weight 
of inverse volume, they become essential when a con- 
densed phase appears since then the occupation of one 
particular channel is proportional to the volume itself. 
The advantage of eliminating only the contributions of 



single channels as proposed in Refs. [l| and [H is that 
the formation of pairs and their condensation can be de- 
scribed within the same approximation. We note that 
the proposed elimination shemes yield similar results as 
the often used mean field approximations with anomalous 
functions [7] . Further improvements of the T-matrix can 
be achieved if the interference of particle-particle chan- 
nels and particle- hole channels are treated as noted e.g. 
in [||. 

From these experiences which gives us the possibil- 
ity to treat the condensed and non-condensed phase on 
the same footing we expect to gain more insight also to 
the problem of interacting Bose systems. As a promi- 
nent example, currently much debated, we present the 
application to these systems and show the appearance 
of 'inaccessible' condensed phases. The theory predicts 
that a Bose gas with contact interaction and a scattering 
length ao < 0, i.e., attractive interaction, undergoes the 
so called Evans-Rashid transition 0, [l(| , which is the bo- 
son analogue to the BCS transition for fermions. It has 
been shown by Stoof 11011 - however, that such Bose gas 
becomes liquid or solid [iCj, [llj before the Evans-Rashid 
transition or Bose condensation can appear. On the other 
hand a BCS-like molecular Bose-condensation as well as 
a condensation of both atoms and molecules is especially 



predicted near a Feshbach resonance Here 
we find three quantum phases: (I) An ideal gas with or 
without bound states, (II) a gaped phase and (III) the 
Bose-Einstein condensate. Whether the latter two phases 
are experimentally accessible is however unclear so fare. 



Let us assume a homogeneous Bose gas with an at- 
tractive interaction strong enough to form a two-particle 
bound state. Such bound states can be found as poles of 
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the T-matrix, which reads in ladder approximation 
T^wj»,fc)=^,fc)-i53v g (p J fc')a,(w,* , )T 9 (w ) * / ,*),(l) 



i.e., the Bethe-Salpeter equation, where q is the center of 
mass momentum and 

G q (iz,p) = T'^G (iw,p)G (iz - iuj,q - p) (2) 



is the two-particle propagator. The Matsubara frequen- 
cies u> are even multiples of 7rT and fl is the system vol- 
ume. By closing the loop, including exchange and sum- 
ming over all q we obtain the selfenergy 



S(iw, k) 



T 



ft X! [ 7 "?(* Z ' k > fc ) + T S Z i fc ' 1 ~ 



z.q 



xG(iz — iu), q — k). 



(3) 



For fermions this is called the Galitskii-Fcynman approx- 
imation. We use the sign convention of Ref. [lq . The 
set of equations is closed by the Dyson equation 



G = G° + G°EG 



(4) 



with the free Green function, G° = — eo(k)] and the 

free single-particle energy eo(k) = H 2 k 2 /2mo — Mo- 
Let us assume a bound state with momentum q and 
energy oj q , i.e., l/T q (uj q — 2{i,p, k) = 0. The chemical po- 
tential /i will be specified later. We further assume that 
the system has excitations at energies E^ , represented by 
poles of the dressed Green function, i.e., 1/G(E^, k) = 0. 
With the Feynman trick we convert the sum over Mat- 
subara frequencies in Eq. ([3]) into a sum over the poles of 
the T-matrix and the dressed Green function [171 ] which 
leads to 

E(iu,k) = J2A 2 q (k)G{uj q -2ti-iuj,q-k) 



where 



+T 9 (£ g ± _^w,fc, ( 7-fc) 



f B (Uq - 2/i) 
ft 



(5) 



xRes [T q (z, k, k) + T q (z, k,q - k), uj q — 2/i] (6) 



and 



4 = ^^^{G(z,k),E±] 



(7) 



The Bose distribution is /s(e) = ]exp(e/T) — l] -1 , the 
residuum of function h(z) at its pole e is denoted by 
Res[/i(z) , e] . The chemical potential /i increases with den- 
sity until it approaches the lowest bounded state, which 



is the state of zero momentum, q — 0. At this point, 
2/i = u>o, the Bose distribution becomes divergent with a 
value proportional to the system volume, i.e., the Evans- 
Rashid transition which means the onset of a conden- 
sation of zero momentum bound states. The divergent 
Bose distribution results in a singular contribution of the 
(q, ui q — 2/i) = (0, 0) channel to the selfenergy 



E s in(iw, k) = A (k)G(—iu, —k). 



(8) 



If the interaction allows for several bound states 
the singular contribution will always come from the 
one with lowest energy. All the other contribu- 
tions to the selfenergy are regular being propor- 
tional to We use a quasi-particle approximation 
Ti leg (iuj, k) = Y,(iw, k) — E s j n (za;, k) « E reg [e(&:), k] with 
the quasi-particle energy 



e(k) = eo(fc) + £ rog [e(fc), k] 



K 2 k 2 
2m 



V + 0(k 3 ). (9) 



Expansion for small momenta determines an effective 
chemical potential [i and an effective mass m. The quasi- 
particle damping could be considered as well if the imag- 
inary part of the self energy is taken into account [18j |. 
For the present message we neglect this effect. 

As argued in Refs. [lj and [5[, the Galitskii- Feynman 
approximation is not able to give the correct dispersion 
in the condensed phase. To avoid this problem the loop 
of the singular self energy contribution © has to be con- 
structed with reduced Green functions which do not in- 
clude the singular self energy 



G ^ — G° + G°E rcg G^ 



1 



iu; — e(k) ' 



i.e., 



Y, s i n (iLu, k) = A (k)Gfi(— iw, —k) 



(10) 



(11) 



and also the two particle propagator of the (0, 0) channel 
has to be constructed with one reduced line 

So(0,p) = G (»w,p) G^ {-iu, -p) . (12) 

The dressed Green function follows from the Dyson equa- 
tion flU), Eqs. (HUD and (JTTJ as 

G = Ga + GE S inGA ~ ' , t — ^-=5, (13) 
(iu>y - E 2 

with a two-branches quasiparticle dispersion E k = ±Ek , 



E k = 



^Je 2 (k)-A 2 (k) 



(14) 



This presents a closed set of equations again. With the 
Green function (|13|) we can now calculate the total den- 
sity 



i^coth^-I 

2E k 2T 2 



(15) 
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For £o = the ground state will give the dominant con- 
tribution to the total density, i.e., there is a condensation 
of quasiparticles. A comparison with (fbT| and ^ yields 
the condition ^ = [i — £ reg [e(0),0] = — A n (0) which 
corresponds to the Hugenholtz-Pines theorem [19j, the 
condition for the usual Bose condensation. 

So far we have only applied the Galitskii-Feynman 
and a quasi-particle approximation. To calculate loq, 
fi and Ao(fc) further approximations and the specifi- 
cation of the interaction are necessary. For the sake 
of simplicity we assume an s-type separable interaction 
Vg{p, k) = A<?p_|<7fe_ j with Yamaguchi form factors [20| . 
9p = 7 2 /(7 2 + P 2 ): where the range of the potential is 
controlled by 7 while A represents the strength of the in- 
teraction. We also neglect the higher orders in k on the 
right hand side of ©. For a separable interaction the 
T-matrix can be expressed as 



H = -5 e 



T q {iz 7 p,k) 



k< 



L G q (iz, k') 



(16) 



From Eq. © follows A (fc) = Ag k . The T-matrix fl6]) 
diverges giving a bound state only if A < 0, therefore the 
condensation of bound states as mentioned above is pos- 
sible only for attractive interaction. Now we are ready 
to perform the limit of infinite volume f2. The momen- 
tum sums will turn into integrals J2k ~ > ^ J^ 3 ^/(2 7r ) 3 J 
and the singular contributions have to be considered sep- 
arately. The condensate density Uq and quantities ujq, /i, 
A can be determined for a given temperature, total den- 
sity and interaction strength with the help of two sub- 
sidiary functions: 



.(A*, A) = 



AbOo,M, A) = - 



d 3 k 
(2tt) 3 

d 3 k 



^coth^ 
2E k 2T 



(17) 



(18) 



Here n gas is the density of the non-condensed Bose gas 
and At, is the interaction strength which would create a 
bound state of energy ujq. We will use these functions to 
discuss the phase diagram. 

As already mentioned there are three phases (I-III) . 
The single-particle dispersions (TTJJ at phases (I-III) are 
shown in Fig. [T] Critical lines between the low, medium 
and high density phases (I), (II) and (III) are speci- 
fied in the phase diagram of Fig. O Figure [3] shows 
the density dependence of cjo, A* an d A for A = — 5A c o- 
In the following we discuss these results in more de- 
tail. According to Eq. © the A vanishes as I/O for 
ujo > The dispersion (fl4)) is thus quadratic lead- 

ing to the ideal gas behavior. Since Eq = — /i > 
there is no Bose condensate, no = 0. The chemi- 
cal potential is determined from the total density n 
via n = n gas (/i, 0). The binding energy follows from 
A = Ab(wo,/i, 0). The two-particle propagator in this 




FIG. 1: Single-particle dispersion for the three different cases: 
(I) quasi ideal, (II) gap, (III) gap-less and linear disper- 
sion. The energy scale is given by the interaction range, 
e 7 = h 2 -y 2 /2m. 
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FIG. 2: Phase diagram for the three different cases. The 
critical densities are given in units of the ideal critical density 
for Bose condensation nid and the couplings strength in units 
of the critical interaction strength in vacuum A c o. Above n c b 
the formation of bound pairs is possible. At n cg appears the 
transition from the quasi ideal phase (I) to the Evans-Rashid 
phase (II) with condensate of bound pairs. Above n cc there 
is an additional condensation of quasi particles, i.e., the usual 
Bose condensation phase (III). 



case is Qo(iui,k) = coth[e(fc)/2T]/[2e(fc) — ioj]. In vac- 
uum the bound state vanishes for A > — A c n = — 

Lu 7717 

For |A| < A c o the bound state appears at the critical 
density n c b given by A = Afe[0, n{n c h), 0], The criti- 
cal density n c b is the dashed line in the phase diagram 
Fig. As already mentioned lvq and 2/j, move towards 
each other with increasing n. At luq = 2\x the conden- 
sation of bound states starts driving the Evans-Rashid 
transition [^, [l(| • The critical density of the transition 



'■eg 



follows from A = A&[2/z(n cg ), fj,(n cg ), 0]. In the phase 
diagram of Fig. [2] the critical density n cg is denoted by 
the solid line. Above n cg a finite A appears due to the 
condensation of bound pairs. The chemical potential /1 
is then pinned to the pair condensate u>q — 2fi. The 
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FIG. 3: The binding energy, chemical potential and A in the 
three different regions according to a cut through the phase 
diagram in Fig. [2] at A = — 5A c q. 



dispersion now consists of two branches separated by a 
gap 2\J y? — A 2 . The gap implies no = 0. The values 
of u and A follow from the equations n = n gas (/^,A) 
and A = Af,(2/i, a, A). In this case only the two-particle 
propagator ([HD, g o {0,p) = coth(E p /2T)/2E p , is needed. 
With increasing density, A increases and [i decreases. As 
can be seen in Fig. [31 their absolute values approach each 
other, i.e., the gap decreases. At u = —A the gap van- 
ishes and the condensation of quasiparticles starts. The 
critical density of this transition is n cc — n gas (— A, A) 
with A given by A = A b (-2A, -A, A). Due to the 
condensation of quasiparticles the zero momentum con- 
tribution to the density (TPo)) has to be split off and 
has to be considered separately as condensate density 
no- The chemical potential is pinned to the conden- 
sate, fj, = —A and uio = — 2A. In analogy to the den- 
sity also the zero momentum component of the sum 
in (fTB]) has to be split off. Therefore no and A fol- 
low from the equations n = no + n gas (— A, A) and 
1/A = -n / A + 1/A b (-2A, -A, A). The dispersion fH]) 
becomes gap-less and linear for small momenta, 



E k = J — 



2A 2 



m H 2 "f 2 



hk + 0{k 2 



(19) 



In leading order of the interaction, i.e., for small A and A, 
we have A = — Xno- In the limit of contact interaction, 
7 — > oo, one then recovers the well known Bogoliubov 
dispersion 



E k 



fh 2 k< 



V 2m 



-Ann 



H 2 k 2 



-Xn 



m 



hk + 0(k 2 ) (20) 



leading to a positive and real sound velocity. 

Finally we have to check whether the system is stable 
or not. To this end we present in Fig. 0] the pressure 



FIG. 4: The pressure versus inverse density in the three dif- 
ferent phases compared to the ideal gas behavior, p — nT. 



TABLE I: Interaction parameters and the typical energy and 
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-1.057 
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0.95 



evaluated via 



P2~Pi= / n[n',A(fi')]dfi'. 



(21) 



The pressure shows a kink at the Evans-Rashid transition 
which is therefore of second order as well as the transition 
to Bose condensation of atoms. Due to the attractive in- 
teraction the pressure decreases above the Evans-Rashid 
transition density n cg . The interaction potential does not 
include short range repulsion or a hard core which leads 
to the artifact that the decrease of the pressure cannot be 
compensated. As the density increases further the result- 
ing pressure becomes negative leading to the collapse of 
the system. Assuming a Van-der-Waals behavior of the 
pressure one would expect the condensation into a liquid 
to start already below n cg . The exact description of this 
condensation line is still open, however. Naturally, the 
theory of three phases described above is plausible only 
for positive pressure. On the other hand the usual densi- 
ties in the experiment are so low, that the condensation 
into a liquid can be avoided during the experiment. 

We restricted ourselves in this paper to the case of 
bosons with zero total spin, e.g. 1 H and 4 He. To ob- 
tain the interaction parameters A and 7 as shown in Ta- 
ble U a comparison with numerical or experimental data 
is necessary. For 4 He A and 7 were derived form the 
scattering length an = 93A [21( and the effective range 
ro = 7.298A 21 1 . In this way one can also obtain the 
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vacuum binding-energy woo = 129.6neV and bond length 
(r) = 49A, which are in good agreement with the nu- 
merical and experimental data w o = (130 ± 7)neV (2lj 
and (r) = (52 ± 4)A respectively. For 1 H the dis- 
sociation energy Dq = 4.478eV [23| and the bond length 
(r) = 0.7414A [23 were used to fit the potential describ- 
ing the 1 H 2 ground state. 

Our results are very similar to those obtained by Stoof. 
We find two phase transitions. The first one is the Evans- 
Rashid transition between a quasi-ideal gas and a BCS- 
like phase connected with a condensation of bound states. 
The second phase transitions is the onset of the usual 
Bose condensation of quasiparticles. We obtain the ex- 
pected linear dispersion in the Bose condensate region 
without the use of anomalous propagators. The kinks in 
the chemical potential and the pressure indicate that the 
Evans-Rashid transition as well as the transition to the 
phase of condensed atoms are of second order. The de- 
crease of the pressure above n cg and the change of its sign 
are a strong signatures of a gas-liquid or gas-solid transi- 
tion and therefore the quantum phase transitions might 
be inaccessible. Concluding we have employed a novel 
many-body correction scheme for the T-matrix which al- 
lows us to describe multiple condensed phases as well 
as gaped phases simultaneously on the same theoretical 
foundation. 
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